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Objectives 


To  develop  efficient  and  robust  ways  to  match  solutions  of  the  Reynolds- Averaged  Navier-Stokes  equations 
to  the  large-eddy  simulation.  These  methods  will  be  tested  in  flows  with  pressure-driven  and  unsteady 
separation  and  in  three-dimensional  boundary  layers. 


Summary  of  results 


We  have  demonstrated  the  applicability  of  synthetic  turbulence  coupled  with  the  cont rolled-forcing 
method  to  generate  eddies  at  the  interface  between  RANS  and  LES  in  hybrid  calculations.  By  investi¬ 
gating  the  effects  of  the  model  parameter,  and  developing  guidelines  to  set  their  values,  we  obtained  shorter 
transition  regions  and  improved  model  accuracy.  It  was  found  that  realistic  turbulence  (with  the  correct 
statistics)  can  be  generated  within  5  bound arv- layer  thicknesses  of  the  RANS/LES  interface  even  in  cases 
(such  as  the  accelerating  boundary  layer)  in  which  the  data  supplied  by  the  RANS  is  inaccurate,  or  the 
assumptions  on  which  the  forcing  is  based  are  invalid. 

We  applied  this  technique  in  a  variety  of  turbulent  Hows  including  phenomena  such  as  freestream  accel¬ 
eration,  separation,  and  mean  How  three-dimensionality,  with  uniformly  good  results.  A  single-block  hybrid 
calculation  was  then  performed  in  a  geometry  that,  although  simple,  presented  several  difficulties:  because  of 
the  shallow  pressure-driven  separation,  incorrect  prediction  of  the  upstream  flow  results  in  significant,  errors 
that  propagate  downstream,  as  the  shear  layer  instability  acts  as  an  error  amplifier.  We  demonstrated  that 
unless  turbulent  eddies  arc  artificially  generated  at  the  RANS/LES  interface,  very  significant  errors  appear 
in  the  How  statistics  even  at  low  order  (skin-friction  coefficient,  mean  velocity  profile).  Establishing  realistic 
turbulent  eddies  capable  to  transport  momentum  energy  and  mass,  appears  to  be  a  critical  factor  for  the 
accurate  prediction  of  shallow  separation  by  hybrid  RANS/LES  methods. 


20090429203 


1 


Accomplishments 


1  Introduction 

1.1  Motivation 

Recent  years  have  seen  the  increased  availability  of  inexpensive  parallel  computer  clusters  able  to  reduce 
the  computational  time  for  the  solution  of  numerical  models.  This  development  has  allowed  to  extend  the 
numerical  solution  of  the  Navier-Stokes  equations  to  turbulent  flows  in  more  realistic  configurations  than 
possible  until  now. 

Turbulent  flows  can  be  studied  by  direct  numerical  simulations  (DNS)  of  the  Navier  Stokes  equation, 
in  which  the  grid  resolution  is  sufficient  to  resolve  all  the  scales  of  motion.  This  method  constitutes  the 
conceptually  simplest  approach  to  the  problem  of  turbulence.  Practically,  however,  the  cost  of  DNS  confines 
this  method  to  low- Reynolds- number  applications:  Assuming  that  computer  power  will  increase  by  a  factor 
of  five  every  five  years,  Spalart  [36]  estimated  that  DNS  will  not  be  applicable  to  the  flow  over  an  airliner  or 
a  car  until  2080. 

In  large-eddy  simulations  (LES)  the  target  is  to  simulate  only  the  eddies  containing  the  bulk  of  t  he  energy 
of  the  flow,  while  modeling  the  dissipative  scales.  These  scales  are  generally  anisotropic  and  dependent  on 
the  boundary  conditions  (therefore,  extremely  difficult  to  model);  the  fact  that  the  small  eddies  are  modeled 
reduces  the  cost  of  LES  considerably  compared  with  DNS,  especially  in  free  shear  flows.  However,  when  LES 
are  applied  to  wall- bounded  flows  at  high  Reynolds  numbers,  the  number  of  grid  points  required  to  resolve  the 
energy-carrying  eddies  is  proportional  to  Re1*8,  which  again  makes  this  technique  only  suitable  for  moderate 
Reynolds- number  applications.  Wall-Modeled  LES  (WMLES),  in  which  the  inner  layer  is  modeled  either  by 
the  Reynolds- Averaged  Navier-Stokes  (HANS)  equations  or  through  approximate  boundary  conditions,  may 
alleviate  this  problem  (see  the  review  in  [22]).  Spalart  [36],  however,  estimates  that  even  the  application 
of  WMLES  to  external  flow  of  aeronautical  interest  is  several  decades  away  for  external  aerodynamics  (in 
environmental  or  oceanographic  flows,  in  which  the  Reynolds  numbers  are  comparably  high,  but  the  boundary 
layers  are  thicker,  WMLES  is  already  applicable  to  practical  problems). 

Traditionally,  high- Reynolds- number  flows  have  been  predicted  using  the  Reynolds- Averaged  Navier- 
Stokes  (RANS)  equations.  RANS  model  are  designed  to  be  accurate  in  a  variety  of  design  flows  (such  as 
thin  shear  layers)  but  lose  accuracy  in  the  presence  of  fluid-dynamical  non-equilibrium  (transition  or  re- 
laminarization,  strong  streamline  curvature,  strong  flow  acceleration  or  deceleration).  A  possible  avenue  to 
solve  practical  engineering  flows,  then,  is  by  using  a  combination  of  RANS  and  LES,  in  which  the  RANS 
equations  are  solved  in  quasi-equilibrium  regions  or  where  they  can  be  accurately  tuned,  while  LES  or  WM¬ 
LES  are  used  in  non-equilibrium  regions,  where  the  RANS  models  are  expected  to  fail.  Hybrid  RANS/LES 
techniques  have  emerged  as  a  tool  that  allows  the  simulation  of  complex  fluid  flows  within  reasonable  amounts 
of  CPU  time. 

One  issue  that  arises  when  hybrid  RANS/LES  methods  are  used  is  the  behavior  of  the  flow  in  the 
transition  zone  between  the  RANS  and  LES  regions  and  how  that  region  affects  the  results  downstream.  In 
the  RANS  zone  the  flow  solution  is  either  steady,  or  only  contains  information  on  the  largest  scales  of  motion 
if  mean- flow  unsteadiness  is  present;  most  (or  all)  of  the  Reynolds  shear  stress  is  provided  by  the  turbulence 
model.  In  the  LES  region  on  the  other  hand,  the  resolved  scales  must  supply  most  of  the  Reynolds  shear 
stress,  and  eddies  must  be  present  to  provide  it.  Typically,  a  transition  zone  exists  in  which  resolved  eddies 
are  generated  gradually  and  grow;  in  this  region  the  flow  may  be  un physical.  Different  techniques  to  trigger 
the  instabilities  have  been  developed. 

One  of  the  most  popular  hybrid  RANS/LES  techniques  is  the  Detached  Eddy  Simulation  [38]  (DES).  In 
this  model  the  switch  between  RANS  and  LES  is  determined  by  the  grid  size:  when  the  mesh  becomes  small 
enough  to  resolve  the  energy-carrying  eddies  the  eddy  viscosity  is  reduced.  Modifications  of  this  model  to 
link  the  RANS/LES  transition  to  the  turbulent  viscosity,  thus  avoiding  premature  switch  from  the  RANS 
to  the  LES  region  have  also  been  proposed  and  show  potential  [35].  In  standard  applications  of  this  method 
(Figure  1),  the  thin  attached  shear  layers  are  modeled  by  RANS,  and  only  away  from  solid  bodies  the 
technique  switches  to  LES,  and  the  RANS/LES  interface  occurs  generally  in  the  separated  shear  layers, 
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Figure  1  Sketch  of  the  interface  region  between  RANS  and  LES  regions  in  a  massively  separated  (low. 


where  the  strong  instability  caused  by  the  inflection  points  in  the  velocity  profiles  is  extremely  beneficial  for 
the  generation  of  energy-carrying  eddies:  any  disturbance  present  in  the  flow  is  quickly  amplified,  resulting 
in  a  very  short  transition  region.  However,  an  important  future  held  of  application  for  hybrid  method  is  in 
aeronautical  applications  such  as  that  sketched  in  Figure  2,  in  which  attached  boundary  layers  are  computed 
using  the  RANS  approach,  while  LES  is  used  in  the  separated  flow  regions.  The  application  of  DES  to  flows 
with  weaker  instability  mechanisms  may  create  problems:  if  the  turbulent  eddies  do  not  grow  sufficiently  fast, 
the  simulations  can  become  inaccurate,  as  the  modeled  shear  stress  decreases  while  the  resolved  one  does  not 
increase  sufficiently  fast.  Such  quantities  as  the  mean  velocity  profile,  skin  friction  coefficient  and  separation 
location  can  be  in  error.  In  fact,  when  DES  is  used  as  a  wall-layer  model  [20],  no  such  instabilities  exist, 
and  a  significant  transition  region  is  observed  between  the  smooth  inner  layer  (in  which  RANS  equations 
are  solved)  and  the  outer  LES  region,  in  which  turbulent  eddies  are  gradually  formed  (see  the  discussion  in 
[22]).  This  transition  region  results  in  a  displacement  of  the  logarithmic  layer  and  errors  in  the  prediction 
of  the  skin-friction  coefficient  in  the  plate. 

Under  a  previous  Grant,  the  Pis  investigated  the  generation  of  inflow  conditions  for  large-eddy  sim¬ 
ulations.  The  methodology  developed  can  be  used  for  hybrid  RANS/LES  models  in  which  two  separate 
calculations  are  carried  out:  with  reference  to  Figure  2,  a  RANS  solution  of  the  entire  domain  is  hrst  per¬ 
formed,  which  supplies  the  boundary  conditions  for  the  LES  region.  The  present  grant  aims  to  optimize  this 
method,  and  use  it  in  actual  hybrid  calculations,  stressing  flows  with  pressure-driven  or  unsteady  separation 
(as  opposed  to  the  massively  separated  flows  to  which  DES  has  been  applied  successfully).  In  these  config¬ 
urations,  the  accurate  prediction  of  the  momentum  transport  due  to  the  resolved  eddies  is  critical;  as  we 
shall  show,  the  absence  of  turbulent  eddies  at  the  RANS/LES  interface  may  lead  to  significant  inaccuracies 
in  the  prediction  of  the  separation  point,  and  of  the  flow  downstream  in  the  recovery  region. 


Figure  2:  Sketch  of  the  interface  region  between  RANS  and  LES  regions  for  an  aerodynamic  application. 
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1.2  Organization  of  this  report 

This  report  is  organized  as  follow:  we  begin  by  summarizing  the  forcing  method  that  had  been  developed 
during  the  previous  grant  cycle.  We  will  then  discuss  the  optimization  of  the  model  parameters.  Finally,  we 
will  present  the  results  of  actual  hybrid  calculations  of  the  How  over  a  contoured  ramp. 

2  Turbulence  generation  techniques 

2.1  Review  of  existing  methods 

We  focus  on  a  class  of  applications,  of  the  type  shown  in  Figure  2,  in  which  RANS  is  used  in  equilibrium 
regions  of  the  How,  while  LES  is  performed  in  a  small  part  of  the  domain,  possibly  including  attached 
boundary  layers  as  well  as  separated  How.  In  this  configuration  the  generation  of  energy-  and  momentum- 
carrying  eddies  can  be  expected  to  be  slow,  as  the  instability  mechanism  in  the  attached  boundary  layer  is 
weaker  than  in  the  separated  shear  layers  that  dominate  the  dynamics  of  massively  separated  flows.  It  is 
the  purpose  of  this  work  to  investigate  methods  to  speed  up  the  development  of  such  turbulent  structures, 
to  obtain  a  realistic  distribution  of  eddies  in  the  shortest  possible  length  for  different  fluid  How  problems. 

Reviews  of  turbulence  generation  methods  for  problems  of  this  type  can  be  found  in  [9,  6],  among 
other  places.  Here,  we  summarize  the  main  issues  involved.  Two  methods  can  be  used  to  perform  hybrid 
calculations:  the  RANS  and  LES  zones  can  be  computed  separately,  with  the  RANS  information  used  to 
assign  the  inflow  boundary  condition  for  the  LES;  or,  a  single  calculation  can  be  carried  out  (as  in  I)ES)  in 
which  the  model  switches  from  RANS  to  LES.  The  simplest  method  to  generate  turbulence  at  an  interface 
between  a  RANS  and  an  LES  zone,  if  separate  calculations  are  used,  is  by  prescribing  a  mean  flow  to  which 
velocity  fluctuations  that  meet  certain  criteria  are  superposed.  If  a  single  calculation  is  performed,  the 
same  result  can  be  achieved  by  including  a  forcing  term  in  the  equations  that  generates  the  desired  random 
perturbations. 

Random  fluctuations  with  given  moments  and  spectra  have  been  used,  for  example,  in  [13,  12].  Since  the 
fluctuations  lack  phase  information,  however,  the  turbulence  levels  decays  rapidly,  and  only  some  distance 
downstream  of  the  inflow  plane  (or  of  the  RANS/LES  interface)  the  turbulent  eddies  are  regenerated.  More 
advanced  methods  that  attempt  to  reproduce  the  eddy  structure  of  a  wall-bounded  flow  at  the  inflow  have 
been  proposed  [32,  26.  11,  1,  3,  16].  They  achieve  improved  results  (compared  with  the  simple  superposition 
of  random  fluctuations)  by  a  judicious  assignment  of  the  turbulence  spectrum  and  by  trying  to  match 
the  flow  anisotropy.  Since  none  of  these  method  contains  realistic  phase  information  between  the  modes, 
an  adjustment  region  downstream  of  the  inflow  is  unavoidable.  In  this  region  the  initial  fluctuations  are 
selectively  amplified  or  dissipated  by  the  flow,  and  realistic  turbulent  eddies  are  generated. 

This  was  shown  by  Keating  et  al.  [9],  who  compared  various  types  of  inflow  conditions  for  LES,  among 
them  the  use  of  random  noise,  the  adapted  database  method  by  Schliiter  [27],  the  synthetic  turbulence 
generation  method  proposed  by  Batten  et  al.  [1],  and  a  technique  based  on  controlled  forcing  proposed  by 
Spille-Kohoff  and  Kaltenbach  [39].  The  latter  method  is  based  on  the  selective  amplification  of  strong  bursts 
downstream  of  an  inflow  supplied  by  a  synthetic  turbulence  generation  method.  A  controller  is  used  at 
several  location  downstream  of  the  inflow  to  determine  the  amplitude  of  a  forcing  term  in  the  wall-normal 
momentum  equation.  This  term  reinforces  the  more  realistic  eddies,  by  requiring  that  a  target  Reynolds 
shear-stress  profile  (obtained  from  the  RANS,  or  from  experiments)  be  achieved.  Keating  et  al.  [9]  found 
that  the  adapted-database  and  the  controlled- forcing  methods  resulted  in  significantly  shorter  development 
lengths  than  any  of  the  other  techniques.  Later,  Keating  et  al.  [6]  applied  the  controlled-foreing  method 
in  a  hybrid  RANS/LES  framework  for  the  simulation  of  boundary  layers  in  favorable  and  adverse  pressure 
gradients,  and  again  found  that  it  produced  physically  realistic  turbulence  in  short  distances.  They  also 
observed  that  the  quality  of  the  RANS  data  used  affects  the  results  significantly;  this  is  true  in  particular  in 
favorable  and  ad  verse- pressure-gradient  boundary  layers. 

In  the  following  we  will  first  describe  the  methodology  used,  then  discuss  the  controller  function,  and  the 
optimization  of  its  parameters.  The  result  of  its  application  will  be  presented  next. 
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Figure  3:  Sketch  of  the  geometric  configuration. 


2.2  Methodology 

2.2.1  Problem  setup 

As  mentioned  above,  hybrid  methods  can  be  developed  both  for  a  single  grid  (as  in  the  DES  applications 
[38]),  or  using  separate  grids  in  the  separate  zones.  This  last  approach  allows,  if  desired,  to  decouple  the 
RANS  calculation  from  the  LES  one  (as  is  the  case  in  the  study  of  the  flow  in  a  turbine  bv  Schliiter  et  ai 
[28]).  For  the  purposes  of  testing  eddy-generation  methods,  the  latter  approach  is  the  most  economical,  since 
the  RANS  is  performed  only  once,  and  its  results  can  be  used  both  to  supply  the  inflow  conditions,  and  also 
to  control  the  forcing.  We  will,  therefore,  follow  this  methodology.  The  tests  described  will  all  be  performed 
using  the  same  strategy  as  [9,  6]  (see  Figure  3):  first  we  perform  a  reference  LES  of  the  entire  domain  of 
interest;  then  we  perforin  a  RANS  calculation  of  the  equilibrium  region,  and  use  the  RANS  statistics  to 
assign  the  inflow  of  the  LES.  Comparisons  between  the  hybrid  RANS/LES  and  the  reference  calculation 
allow  ns  to  evaluate  the  effectiveness  of  the  method. 


2.2.2  Governing  equations  and  numerical  method 


In  LES  the  velocity  field  is  separated  into  a  resolved  (large-scale)  and  a  subgrid  (small-scale)  field,  by  a 
spatial  filtering  operation.  In  the  unsteady  RANS  (URANS)  approach,  on  the  other  hand,  the  Reynolds 
decomposition  is  used  to  separate  the  mean  from  the  fluctuating  part  of  the  velocity.  Despite  this  difference, 
the  governing  equations  for  LES  and  URANS  have  the  same  form 
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In  these  equations,  an  overline  denotes  a  quantity  averaged  over  one  grid  cell  and  time  step,  and  rfJ  = 
UiUj  —  utUj  are  either  the  snbgrid-seale  (SGS)  stresses  or  the  Reynolds  stresses  that  must  be  modeled.  In 
both  regions  an  eddy- viscosity  model  is  used  to  close  the  above  set  of  equations: 
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The  model  used  in  the  LES  region  to  parameterize  the  SGS  stresses  is  the  Lagrangian  dynamic  eddy- viscosity 
model  [5  17].  In  the  RANS  region  the  Spalart-Allmaras  (SA)  one-equation  model  [34]  is  used.  The  model 
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Here  \il\  is  the  magnitude  of  the  vorticity,  and  the  function  fw  is  given  by 
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The  constants  in  the  model  are  q> i  =  0.1355,  a  =  2/3,  Cb 2  =  0.622,  k  =  0.41,  cw\  =  c^/h2  +  (1  +  c.b2)/cr , 
cw2  =  0.3,  ^3  =  2.0  and  cvi  =7.1.  For  the  RANS  application,  the  length  scale  d  is  the  distance  from  the 
wall,  yw. 

The  governing  equations  were  solved  on  a  Cartesian  staggered  grid.  Conservative  second-order  finite 
differences  were  used  for  spatial  discretization.  A  fractional-step  method  [10]  with  a  second-order  implicit 
Crank-Nicolson  time- advancement  for  the  wall-normal  diffusion  term  and  a  third-order  explicit  Runge-Kntta 
method  for  the  remaining  terms  was  used  for  time  integration.  Periodic  boundary  conditions  were  used  in 
the  spamvise  direction,  convective  conditions  [21]  at  the  outflow,  and  the  freestream  condition  [15], 
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at  the  upper  boundary;  here  6*  is  the  boundary  layer  displacement  thickness  and  Uoo  is  the  freestream 
velocity.  We  calculated  dS* /dx  using  a  linear  regression  on  the  $*(#)  distribution.  The  Poisson  equation  is 
solved  by  a  direct  solver  based  on  spanwise  Fast-Fourier  Transforms  followed  by  cyclic  reduction.  The  code 
is  parallelized  using  Message  Passing  Interface  routines. 


2.2.3  Turbulence  generation  method 


The  method  used  to  generate  inflow  turbulence  is  the  cont  rolled-forcing  method,  described  in  detail  in  [9,  6]; 
here  we  summarize  its  main  features,  concentrating  on  the  parameters,  and  in  the  following  sections  we  will 
make  important  modification  and  generalization  to  that  model.  The  turbulence  generation  method  consists 
of  two  parts:  the  creation  of  synthetic  turbulence,  and  the  controlled  forcing.  The  synthetic  turbulence 
generation  method  of  [1]  is  used  to  an  unsteady  velocity  field  at  the  inflow  plane  of  the  LES  region.  An 
intermediate  velocity,  is  first  constructed,  using  a  sum  of  sines  and  cosines  with  random  phases  and 
amplitudes: 


(•Tj.  0  =  \J  51  [p>‘  cos(d"J"  +  u>nt)  +  q"  sin(^x"  +  wnf)] 


(8) 


where 


Xj  =  2nXj/  Lb,  t  =  2ni/rb, 


(9) 


are  spatial  coordinates  normalized  by  the  lengt  h-  and  time-scale  of  the  turbulence.  In  the  above,  Tb  =  IC/e 
and  =  7b  14  are  the  turbulence  time-  and  length-scale,  and  14  =  K}!2  is  the  velocity  scale.  The  random 
frequencies  uon  =  iV(l,l)  are  taken  from  a  normal  distribution  N(y,cr2)  with  mean  fi  =  1  and  variance 
a2  —  1.  The  amplitudes  are  given  by 


where  Q ,  =  N(0, 1),  and 
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are  modified  wavenumbers  obtained  by  multiplying  the  wavenumbers,  d™,  by  the  ratio  of  the  velocity  scale 
Vf,  =  Lb/rt,  to  cn,  given  by 
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The  wave-numbers  d™  =  N(Q, 1/2)  are  chosen  from  a  normal  distribution  with  variance  1/2,  resulting  in  a 
three-dimensional  spectrum  that  behaves  like  d4  exp(— d2).  Although  the  wave-numbers  d-1  are  distributed 
isotropically  in  a  sphere,  dividing  them  by  cn  tends  to  elongate  those  wave-numbers  that  are  most  closely 
aligned  with  the  largest  component  of  the  Reynolds-stress  tensor,  and  contract  those  aligned  with  the  smaller 
ones.  This  results  in  a  more  physically  realistic  spectrum  of  turbulence,  with  eddies  that  (near  the  walls) 
are  more  elongated  in  x,  and  tend  to  be  more  spherical  in  the  channel  center.  The  synthetic  turbulent 
fluctuation  field  is  finally  reconstructed  by  a  tensor  scaling: 


u'i  =  (likVk 


(13) 


where  is  the  Cholesky  decomposition  of  the  Reynolds-stress  tensor.  The  method  requires  turbulence  time- 
and  length-scale,  77,  =  JC/e  and  Lb  =  r^/C1  2  (where  /C  is  the  turbulent  kinetic  energy  and  e  the  dissipation) 
that  must  be  supplied  by  the  RANS  solution.  When  the  SA  model  is  used,  following  [2,  18],  we  estimate  k 
and  e  from 

s/C^IC  =  |(-uV)| ;  c  =  C^/ut.  (14) 

The  cont rolled-forcing  method  consists  of  a  forcing  term  added  to  the  wall-normal  momentum  equat  ion 
that  amplifies  the  velocity  fluctuations  in  that  direction,  thus  enhancing  the  production  term  in  the  shear- 
stress  budget.  The  forcing  amplitude  is  determined  by  a  PI  controller,  which  has  two  components:  a 
proportional  part  and  an  integral  one.  The  two  outputs  are  added  together  to  form  the  actuating  signal  of 
the  system  to  control,  the  Navier-Stokes  equation  system. 

The  error  is  defined  as 

e  =  (u'v')tar9€t  -  (u'v'Y  (15) 

where  (ti'v')iar9et (xQ,  y)  is  the  desired  Reynolds  shear  stress  at  the  control  plane  x  =  xQ  obtained  from  the 
RANS  solution  and  which  is  the  target  of  the  control  procedure;  (u'v'Y (xQ,  is  the  Reynolds  shear  stress, 
averaged  over  some  time  interval.  All  terms  are  evaluated  at  the  control  plane  x0,  and  may  be  function  of 
y,  z  and  t.  The  magnitude  of  the  forcing  is  then  set  to 

fCF(x0,y,  z,t)  =  r(y,  z,  t)  [u(x0,y,z,t)  -  {u)1  (x0,y,  z,t)]  (16) 

where 

r(t)  =  Kpe(t)  +  A'/  J  e(t')dt'  (17) 

is  the  output  of  the  PI  controller.  The  significance  of  the  two  constants,  Kp  and  A'/  and  the  optimal  value 
of  the  length  of  the  time  average  window,  will  be  discussed  later.  The  forcing  so  defined  is  added  to  the 
{/-momentum  equation.  Enhancing  the  v'  fluctuations  through  events  with  large  u'  [the  term  in  square 
brackets  in  (16)]  has  the  effect  of  accelerating  the  production  of  Reynolds  shear  stress. 


3  Analysis  of  the  controller  parameters 

3.1  Preliminary  study  in  a  flat-plate  boundary  layer 

Although  the  controlled-forcing  method  has  shown  some  promise,  its  development  so  far  has  been  heuristic. 
The  first  step  of  our  research  consisted  in  investigating  more  closely  the  link  between  the  model  parameters 
and  the  flow  physics,  to  establish  the  robustness  and  stability  of  the  method  and  to  improve  its  performance. 
The  test  case  used  to  evaluate  the  parameters  in  the  model  was  a  flat-plate  boundary  layer;  then,  the 
results  were  extended  and  generalized  to  more  complex  flows  with  different  characteristic  time-  and  length- 
scales,  such  as  a  fiat-plate  boundary  layer  subjected  to  favorable  and  adverse  pressure  gradients,  and  a 
three-dimensional  boundary  layer. 
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The  controller  described  above  has  three  parameters:  the  constants  Kp  and  Kj  and  the  width  Tave  of  the 
averaging  window  used  to  obtain  (tiV)L  In  past  applications  [39,  9,  6]  no  investigation  was  carried  out  in 
which  they  were  systematically  varied.  In  this  section  we  examine  the  transient  and  steady-state  responses 
of  the  flow  to  these  parameters,  and  relate  them  to  physical  properties  of  the  flow,  as  well  as  to  each  other. 
The  target  is  to  choose  them  to  obtain  statistically  steady-state  Reynolds  stresses  that  match  the  desired 
ones  in  a  short  distance  without  destabilizing  the  Navier-Stokes  equation  system.  While  we  cannot  claim 
that  to  have  found  true  optimal  values  of  these  parameters,  we  will  refer  to  our  choice  as  “optimal91  if  it 
gives  realistic  turbulence  in  the  shortest  distance,  and  with  the  shortest  transient. 

We  test  them  in  a  fiat-plate,  zero- pressure-gradient  (ZPG)  boundary  layer.  First,  we  carried  out  the 
reference  LES  on  a  domain  240 5*  x  25 5*  x  25 5*  in  the  streamwise,  wall-normal  and  spanwise  directions, 
respectively;  here  5*  is  the  displacement  thickness  at  the  inflow.  From  this  calculation,  the  time-averaged 
Reynolds  shear  stress  and  the  turbulence  scales  were  extracted  at  x/6*  =  80.  At  this  section  a  FES  of 
dimension  1205*  x  255*  x  255*  with  the  synthetic  turbulence  generation  and  the  controlled  forcing  at  the 
inflow  began.  In  this  way  it  was  possible  to  analyze  the  effect  of  Kjf  Kp  and  Tave  without  the  errors 
introduced  by  a  RANS  model.  The  control  planes  were  distributed  over  a  length  of  27.55*  downstream  of 
the  interface  in  the  LES  domain,  with  a  spacing  between  planes  of  one  boundary-layer  thickness  to  allow 
the  flow  to  re-adjust  after  the  forcing  (test  calculations  in  which  the  forcing  was  distributed  continuously 
had  also  carried  out,  with  no  significant  differences  [7]).  The  inflow  Reynolds  number  (based  on  freestream 
velocity  UG  and  displacement  thickness)  was  Re*6  —  1000.  The  grid  spacings  in  the  streamwise  and  spanwise 
directions  were  Ax/5*  =  1.25  and  A z/S*  =  0.385,  while  64  points  were  used  in  the  wall-normal  direction 
(with  y+in  «  1.0).  The  recycling  and  rescaling  method  of  Lund  [15]  was  used  as  inflow  boundary  condition 
for  the  reference  LES.  This  resolution  is  sufficient  to  ensure  accurate  prediction  of  the  mean  velocity  and 
Reynolds  stresses  [8]. 

PI  controllers  give  a  robust  performance  over  a  wide  range  of  operating  conditions  and  are  widely  used 
due  to  their  easy  implementation.  They  are  frequently  used  when  the  mathematical  model  of  the  system  to 
control  is  complex  or  unavailable,  in  which  case  the  parameter  tuning  can  be  based  on  semi-empirical  rules. 
In  unknown  non-linear  system,  however  the  primary  target  of  the  parameter  tuning  is  the  stabilization  of 
the  controlled  system,  and,  marginally,  the  steady  state  performance.  In  the  cont rolled-forcing  turbulence 
generation  method,  the  system  to  control  is  the  filtered  NS  equation  system  with  a  dynamic  subgrid  model. 
The  target  is  the  Reynolds  shear  stress  obtained  from  the  RANS  solution.  The  output  of  the  system  is  the 
instantaneous  u'v'  correlation,  averaged  using  an  appropriate  filter. 

We  begin  by  discussing  the  moving-average  (MA)  filter.  In  the  classical  terminology  used  in  the  feedback 
control  theory,  this  block  is  called  “Measuring  device”:  the  instantaneous  velocity  fluctuations  u[  are  the 
inputs  to  the  MA  block,  which  outputs  an  appropriate  time-averaged  Reynolds  stress  that  can  be  compared 
to  the  target  obtained  from  the  RANS.  We  use  an  exponentially  weighted  moving-average  filter,  which  places 
more  emphasis  on  the  most  recent  data  available.  At  time-step  k  we  use  the  following  expression: 

(xky  =  (l-  A-)  (18) 

At  is  the  time-step,  and  x  is  the  data  to  average. 

The  value  of  A t/Tave  determines  the  memory  of  the  filter.  We  define  an  attenuation  time  Tj  as  the  time 
after  which  the  contribution  of  the  signal  at  some  time  to  the  averaged  data  is  negligible  (say,  5%  of  the 
original  contribution).  It  is  easy  to  show  that,  in  this  case, 

=  log  0-05  n 

Tave  Tavc  Jog 

since  At  >>  Tave,  (19)  gives  Tave  ~  7^/3.  We  expect  that  Td  should  be  of  the  order  of  a  large-eddy  turnover 
time  (LETOT),  which  results  in  Tave  ~  105*/r/r  in  the  present  calculation.  Values  much  larger  than  this 
result  in  long  transients,  as  the  error  affects  the  input  to  the  forcing  for  a  long  time,  and  the  forcing  does  not 
adjust  rapidly  enough  to  the  present  flow  conditions,  but  is  strongly  affected  by  past  events.  This  is  shown 
in  Figure  4(a)  where  the  domain-averaged  TKE  is  plotted  as  a  function  of  time;  in  all  the  calculations  we 
used  values  of  Kp  and  K[  found  to  give  stability  to  the  whole  system  (see  the  next  subsection),  except  for 
one  case  in  which  we  employed  the  values  used  in  [9,  6]. 
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Figure  4:  (a)  Domain-averaged  turbulent  kinetic  energy  for  Tave  =  l,Kj  =  5,  Kp  =  30: - Tavc  = 

10,  Ki  =  5,  KP  =  30;  ---  Tave  =  100,  Kj  =  1,  KP  =  1;  —  Tave  =  100,  Kj  =  5,  I<P  =  30.  (b)  Cf  for  the 
reference  LES  (o  )  and  the  cases  shown  in  (a). 

A  longer  transient  can  be  observed  with  TaVe  =  100  compared  with  Tave  =  10.  With  the  lower  values 
of  KP  and  Kj  used  in  early  applications  of  this  method  the  transient  is  quite  long,  the  magnitude  of  the 
oscillations  is  significant-  and  How  statistics  require  very  long  averaging  times.  Note  that  the  value  of  Kj 
and  KP  used  in  [9,  6]  are  at  the  limit  of  the  stability.  A  stable  choice  of  Kp  and  Kj ,  on  the  other  hand, 
reduces  the  amplitude  of  the  fluctuations  even  if  a  high  value  of  Tave  is  used.  A  low  value  of  Tave  results  in 
incorrect  prediction  of  the  How:  the  short  averaging  gives  values  of  {u'v'Y  that  are  too  far  from  a  stationary 
sample,  and  may  result  in  excessive  forcing  being  applied  even  if  the  How  has  reached  a  realistic  state. 

In  Figure  4(b)  the  skin-friction  coefficient  Cf  —  2rw/ pU^  (where  rw  is  the  wall  stress)  is  shown  for  the 
same  cases.  Excessively  small  values  of  Tave  result  in  incorrect  steady  state  results,  while  the  other  values 
of  Tave  are  in  reasonable  agreement  with  each  other.  The  case  with  Tave  =  10  (f.e.,  Td  of  the  order  of  one 
LETOT)  gives  somewhat  more  accurate  statistics  and  a  much  shorter  transient  (and  reduced  CPU  costs). 
Unless  explicitly  specified  Tave  =  10  for  all  following  cases. 

On  a  related  note,  it  should  be  mentioned  that  previous  applications  of  the  controlled  forcing  method 
[39.  9,  6]  used  spanwise  averaging  (in  addition  to  the  time-averaging)  to  supply  a  smoother  signal  to  the 
forcing.  The  results  shown  above,  however,  indicate  that  this  requirement  must  be  balanced  by  the  need 
to  have  an  error  that  reflect  more  closely  the  local  and  instantaneous  conditions  of  the  How.  We  performed 
calculations  in  which  no  spanwise  averaging  was  performed  and  compared  them  to  a  simulation  in  which  the 
u'v'  correlation  is  averaged  in  the  spanwise  direction  as  well  as  in  time  (Figure  3.1).  We  observe  more  rapid 
adjustment  of  the  How  towards  the  reference  LES,  as  the  forcing  is  more  responsive  to  the  local  state  of  the 
How. 

We  now  turn  our  attention  to  the  other  parameters  of  the  PI  controller  to  evaluate  their  influence  on  the 
transient  time,  steady-state  results  and  stability  of  the  system.  In  general  a  proportional  controller  affects 
the  rise-time  and  the  transient  time  of  the  dynamic  system:  it  modifies  the  bandwidth  of  the  frequency 
response  of  the  closed  loop  system  and  the  gain  of  the  zero  frequency.  If  Kp  is  large,  the  output  of  the 
integral  block  is  more  sensitive  to  the  high-frequency  variation  of  the  error,  and  the  stability  of  the  system 
can  be  affected.  The  integral  control,  on  the  other  hand,  gives  a  large  gain  at  low  frequency  and  reduces 
the  break  frequency  (the  point  where  the  amplitude  spectrum  of  the  transfer  function  is  zero);  its  effect  is 
to  eliminate  the  steady-state  error  and  attenuate  the  high-frequency  disturbances.  In  fact,  the  output  of  the 
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Figure  5:  Skin  friction  coefficient  C/  with  Kj  =5  and  Kp  =  30: - (uV)*  is  spanwise- averaged; 

{u'v'Y  is  not  span  wise- averaged. 
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Figure  6:  Domain-averaged  turbulent  kinetic  energy,  (a)  Kj  =  5  and - Kp  —  0;  15; - 30; - 

500;  (b)  Kv  =  30  and -  Kt  =  1.  ---  5,  —  20,  30. 

integral  control  will  change  over  time  as  long  as  an  error  exists,  but  the  phase  lag  between  the  input  and 
output  of  the  integral  block  increases  for  all  of  the  frequencies  (the  phase  of  the  integrator  block  starts  at 
—90°).  If  Ki  is  large,  the  phase  lag  is  extended  to  higher  frequency,  so  the  system  can  became  oscillatory 
and  potentially  unstable. 

Figure  6  shows  how  the  transient  is  affected  by  Kp  and  Kj.  If  Kp  =  0  (only  the  integral  control  is 
activated)  the  calculation  diverges:  in  the  initial  transient  the  forcing  is  proportional  to  e(t')dt'  and  the 
initial  time-step  is  small  because  of  the  initial  unphysical  flow  in  the  computational  domain.  Because  of 
the  small  time-step  the  integral  action  is  weakened,  so  the  integral  control  needs  a  longer  time  before  it 
can  give  a  correct  control  signal.  As  the  simulation  advances,  however,  the  time-step  is  further  decreased 
because  no  fast  correction  is  present  in  the  domain,  so  the  correct  integral  action  is  even  more  delayed  until 
the  simulation  goes  unstable.  The  case  of  high  Kp  is  the  opposite:  as  soon  as  the  simulation  starts  there 
is  the  fast  correction  due  to  the  proportional  control  Apart  from  the  length  of  the  transient  however,  the 
steady-state  results  did  not  differ  much  in  all  the  converged  calculations,  showing  little  sensitivity  of  the  flow 
to  the  proportional  controller  (at  least  for  Kp  >  1). 

Similarly,  we  found  instability  for  Kj  >  30,  as  can  be  expected  by  the  theory  of  the  PI  feedback  control. 
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Figure  7:  (a)  Instantaneous  and  (b)  integrated  error  for  I\p  =  30  and  Ki  —  5  (lines)  and  Kj  =  20  (lines 

with  symbols). - y+  =  13;  -  y/5 *  =  4.5.  The  curves  for  Kj  —  20  are  shifted  upwards  bv  0.005  in 

(a)  and  0.1  in  (b). 


If  K i  was  close  to  zero,  only  the  proportional  control  was  activated  and  the  steady  state  error  (the  difference 
between  the  desired  Reynolds  stress  and  the  one  obtained  from  the  calculation)  was  large.  Figure  7  shows 
the  error  and  its  integral  at  two  locations,  one  in  the  wall  layer,  the  other  in  the  outer  region  of  the  boundary 
layer,  for  Kp  =  30  and  two  values  of  Kj.  Increasing  Kj  reduces  the  amplitude  of  the  dominant  frequency  of 
the  error,  as  well  as  its  amplitude,  at  least  in  the  outer  layer.  Note,  however,  that  the  forcing  signal  (which 
has  the  integral  error  multiplied  by  A'/)  does  not  change  its  magnitude. 

As  a  result  of  the  tests  described  in  this  section,  we  conclude  that  some  localization  of  the  error,  both 
in  time  and  space,  is  desirable.  In  fact,  inclusion  of  non-local  (in  space  or  time)  effects  may  potentially  give 
instability  to  the  entire  system.  In  that  sense,  the  use  of  distributed  forcing  as  in  [7]  may  result  in  longer 
development  lengths  because  it  uses  non  local  information.  Compared  with  previous  work,  we  obtained 
improved  results  using  a  time-averaging  window  with  a  time-scale  matching  the  LETOT  of  the  flow,  and 
removing  the  spanwise  averaging.  We  observed  that  the  controller  coefficients  mostly  affect  the  length  of  the 
transient  and  the  stability  of  the  system;  for  a  wide  range  of  both  Kj  and  Kp ,  however,  the  flow  statistics 
are  fairly  insensitive  to  the  parameter  values.  In  the  next  section  we  will  apply  the  method  with  the  time- 
averaged  window  using  the  LETOT  as  a  time-scale  to  different  flows,  to  determine  its  performance  (with 
the  new  coefficients)  in  actual  cases. 

3.2  Applications 

3.2.1  Adverse-pressure-gradient  boundary  layer 

Hybrid  simulations  were  next  carried  out  in  an  ad  verse- pressure-gradient  (APG)  boundary-layer  that  un¬ 
dergoes  separation.  The  configuration  of  these  simulations  is  similar  to  that  of  [6],  with  an  inflow  Reynolds 
number,  Res *  =  1260,  and  a  profile  of  imposed  at  the  top  boundary  that  results  in  a  significant  deceler¬ 
ation  of  the  flow  (see  Figure  8).  Due  to  the  strong  adverse  pressure-gradient,  the  flow  separates;  a  favorable 
pressure  gradient  then  closes  the  recirculation  bubble.  The  Reynolds  number  of  this  calculation  is  lower  than 
that  of  the  DNS  by  Na  and  Moin  [19],  but  calculations  carried  out  at  the  same  Reynolds  number  showed 
excellent  agreement  with  the  DNS  data. 

The  computational  domain  used  in  the  full-domain  LES  was  380 5*  x  64£*  x  20 5*  (a  subscript  o  refers  to 
values  evaluated  at  the  inflow  location).  This  LES,  which  used  the  rescaling/recycling  method  at  the  inflow. 
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Figure  8:  Freest  ream  velocity  for  the  adverse- pressure-gradient  calculation. 


had  a  grid  of  384  x  192  x  64.  The  RANS  calculations  for  the  hybrid  cases  were  extended  to  cover  the  entire 
How  domain  t  o  remove  difficulties  of  placing  outflow  boundary  conditions  close  to  reversed-flow  regions;  both 
SA  and  1C  —  e  models  were  used  for  the  RANS  solution.  The  LES  region  started  at  x/S *  ^  80,  and  used 
the  RANS  data  at  that  location  for  the  synthetic  turbulence  generation,  as  well  as  target  for  the  controlled 
forcing  algorithm.  The  S*/uT  obtained  from  the  1C  —  e  model  at  the  interface  location  was  approximately 
110,  which  agrees  well  with  the  value  from  the  full  LES.  To  match  this  value,  Tave  was  fixed  to  38,  so  that 
the  memory  Td  was  approximately  113.  Kj  was  fixed  at  5  and  Kp  to  30. 

Figure  9  shows  the  C/,  mean  velocity  profiles  and  Reynolds  shear  stresses  at  three  locations.  As  the 
boundary  layer  is  subjected  to  the  adverse  pressure-gradient,  C/  decreases  until  the  How  separates  at  x/S*  % 
170.  A  weak  separation  bubble,  which  has  a  length  of  approximately  80(5*,  can  be  observed  in  the  full-domain 
LES.  The  results  from  the  RANS  simulation  in  terms  of  Cf  are  close  to  those  of  the  full  LES.  Switching 
to  the  LES  is  however,  beneficial  (the  prediction  of  Cf  and  the  dimensions  of  the  separation  bubble  are 
both  more  accurate).  Shortly  after  separation  (the  second  profile  in  Figure  9)  some  differences  between  the 
hybrid  cases  can  be  observed,  whereas  inside  the  separation  bubble  the  LES  gives  similar  results.  In  this 
flow,  the  amplification  of  turbulence  in  the  separated  shear  layer  acts  as  a  powerful  mechanism  to  accelerate 
the  generation  of  realistic  eddies. 

Keating  et  al.  [6]  observed  a  much  stronger  dependence  of  the  hybrid  results  on  the  accuracy  of  the 
RANS  model  (see  Figure  11  in  [6]).  This  effect  is  due  to  the  use  of  suboptimal  PI  controller  parameters  and 
averaging  window.  With  the  more  localized  averaging  used  here  the  initial  decay  of  the  wall  stress  that  was 
observed  by  Keating  et  al.  [6]  does  not  occur. 

Figure  10  shows  contours  of  streamwise  velocity  fluctuations  in  a  plane  near  the  wall  for  the  reference 
LES  and  a  hybrid  case.  We  observe  a  very  rapid  development  of  a  physical  streaky  structure.  The  onset  of 
separation  (which  is  strongly  affected  bv  the  flow  state  immediately  before  the  separation  occurs)  is  predicted 
well,  and  the  length-scales  of  the  flow  in  the  separated- flow  region  are  remarkably  similar,  confirming  again 
the  robustness  and  effectiveness  of  the  method. 

3.2.2  Favorable  pressure-gradient  boundary  layer. 

Simulations  of  a  boundary  layer  subjected  to  a  favorable  pressure-gradient  (FPG)  were  performed  next.  In 
this  flow,  if  the  acceleration  is  rapid  enough,  re-laminarization  and  re-transition  of  the  flow  may  occur.  The 
case  studied  here  matches  the  experiments  of  [41].  The  freestream  acceleration  (from  U ^  =  1  to  U0 0  ~  3) 
begins  at  approximately  100(5*  and  is  completed  by  450(5*.  The  flow  acceleration  is  achieved  by  imposing 
a  streamwise  velocity  profile  U^x)  at  t  he  top  boundary  of  the  domain  [15]  Its  magnitude  was  calculated 
from  the  acceleration  parameter,  K  =  {v/U^^dJJ /dx),  experimentally  obtained  from  [41].  The  acceleration 
parameter.  A",  and  the  free-stream  velocity  distribution  Uoo(x)  are  the  same  as  [6].  Since  K  exceeds  the 
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Figure  ():  Adverse- pressure-gradient  boundary  layer,  (a)  Skin  friction  coefficient  Cf  (b)  mean  velocity 
profiles  and  (c)  |>roHles  of  the  Reynolds  shear  stress  at  the  locations  indicated  bv  a  vertical  line  in  part  (a). 
o  Reference  LES; - SA  HANS: - Hybrid  SA  RANS/LES. 


Figure  10:  Streamwise  velocity  fluctuations  in  the  APG  boundary  layer,  in  the  y/6*  =  0.09  plane,  (a) 
Reference  LES;  (b)  hybrid  JC  —  r/LES;  the  red  rectangle  indicates  the  region  where  the  control  is  active. 
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Figure  11:  Favorable-pressure-gradient  boundary  layer  (a)  Skin  friction  coefficient  Cj  (b)  mean  velocity 
profiles  and  (c)  profiles  of  the  Reynolds  shear  stress  at  the  locations  indicated  by  a  vertical  line  in  part  (a). 

o  Reference  LES; .  SA  RANS;  -  hybrid  SA/LES,  Tav€  from  LES  data; - hybrid  SA/LES,  Tave 

from  RANS  data;  . 


critical  value  for  re-laminarization  (Kent  —  3.5  x  10  6)  for  an  extended  region  of  the  flow,  turbulence  is 
expected  to  be  damped  in  the  region  of  high  acceleration.  The  How  then  re-transitions  once  the  acceleration 
is  removed. 

A  full-domain  LES  was  first  performed  using  the  rescaling/ recycling  method  at  an  inflow  Reynolds 
number,  Res*  =  1260.  For  this  simulation  the  domain  length  was  4765*,  the  width  was  205*  and  the  height 
was  205*.  A  somewhat  coarse  grid  (512  x  64  x  64)  was  used,  since  the  comparisons  were  primarily  being  made 
between  this  full  LES  and  a  hybrid  RANS/LES  simulations.  A  previous  LES  study  using  a  finer  grid  has 
shown  excellent  agreement  with  the  experiments  for  this  How  [4];  on  the  present  grid,  the  qualitative  behavior 
of  the  How  was  captured  accurately,  although  the  re-transition  following  the  release  of  the  acceleration  was 
slightly  delayed.  One  hybrid  calculations  was  performed  that  included  a  RANS  domain  3505*  long,  and  an 
LES  region  that  started  at  x/S*  =  225  and  extended  to  4765*  (i.e.,  had  a  length  2515*).  Synthetic  turbulence 
with  controlled  forcing  was  introduced  at  the  RANS/LES  interface. 

The  SA  RANS  model  gives  significant  modeling  error,  being  unable  to  predict  the  relaminarization  and 
the  re-transition  of  the  accelerating  boundary  layer  correctly.  Furthermore,  the  integral  time-scale  at  the 
interface  location  was  approximately  2  5*/uTi0,  while  the  value  obtained  from  the  full  LES  calculation  is 
twice  as  large  in  the  region  where  the  control  was  applied.  Two  hybrid  calculation  were  performed,  one 
that  used  a  value  of  Tt 'd  close  to  the  time-scale  predicted  by  the  RANS  (giving  TaVe  =  1),  and  another  that 
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Figure  12:  3D  boundary  layer.  Flow  angle  a  —  tan  1  U/W.  o  Reference  LES. .  SA-RANS; - hybrid 

RANS/LES,  error  based  on  (uV); - hybrid  RANS/LES,  production- based  error. 


matched  the  LES  time-scale  ( Tave  ~  2.4).  Figure  11  compares  the  results  of  the  various  simulations.  The 
mean  data  at  the  interface,  which  is  supplied  by  the  RANS,  is  not  accurate;  the  controller  tries  to  drive  the 
LES  towards  the  RANS  region,  but  when  Tave  is  small  a  smooth  (uV)  cannot  be  obtained,  and  excessive 
forcing  is  applied,  resulting  in  an  initial  overshoot  of  the  skin- friction  coefficient .  Once  the  control  is  released, 
the  flow  adjusts  fairly  rapidly  towards  the  full-domain  LES.  When  a  longer  average  is  used,  the  controller  is 
more  successful  in  driving  the  LES  shear  stress  towards  the  desired  distribution  (Figure  11(c))  and  the  LES 
adjusts  more  rapidly  and  with  no  overshoot.  Note  that  this  is  a  very  challenging  test  case:  RANS  models 
cannot  be  expected  to  predict  relaminarization  accurately;  furthermore,  since  the  energy  of  the  fluctuations 
is  decreasing  because  of  the  acceleration,  adding  energy  through  the  forcing  term  is  moving  the  flow  away 
from  its  equilibrium.  Despite  tin'  fact  that  this  can  be  considered  an  off-design  application  for  the  controlled 
forcing,  the  results  are  altogether  acceptable. 


3.2.3  3-D  Boundary  layer 

Simulations  were  then  performed  on  a  3D  boundary  layer  obtained  by  applying  a  spanwise  pressure  gradient 
to  a  flat-plate  boundary  layer.  The  pressure  gradient  increased  from  zero,  at  a?/J*  =  100,  up  to  1.5  x  1()~3 
at  a:/(S*  =  150,  and  remained  constant  thereafter.  The  magnitude  of  the  pressure  gradient  was  set  in  such  a 
way  as  to  turn  the  flow  by  45°  near  the  wall  and  24°  at  the  free  stream  bv  the  end  of  the  domain.  Figure  12 
shows  the  turning  angle  near  the  wall  and  in  the  freestream. 

Once  again,  two  simulations  were  carried  out:  a  full-domain  LES,  and  one  hybrid  RANS/LES  calculations. 
The  computational  domain  of  the  reference  LES  was  29 1£*  x  206*  x  20<S*  with  a  grid  of  720  x  100  x  292.  An 
inflow  plane  from  a  precursor  simulation  was  used  with  an  inflow  Reynolds  number  of  Res •  =  1260.  Because 
of  the  turning  of  the  How,  the  grid  had  to  be  refined  in  the  streamwise  direction  in  the  downstream  area.  We 
define  a  local  coordinate  frame  with  £  in  the  flow  direction  and  £  in  the  lateral  direction  in  the  wall  plane: 

A£+  =  Ax+  cos  a  +  As+  sin  a;  A£+  =  sin  a  4-  Az+  cos  a.  (20) 


We  maintained  streamwise  and  spanwise  grid  resolutions  of  A£+  <  50  and  A£+  <  13.  Here,  a  is  the  flow 
angle  at  the  wall  and  wall  units  are  defined  using  the  resultant  wall  stress 
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In  the  hybrid  cases,  the  RANS  equations  were  solved  using  the  Spalart-Allmaras  model.  The  LES  region 
started  at  oo/S*  «  200  where  the  flow  angle  was  32°  near  the  wall  and  12°  in  the  free-stream.  The  SA  model 
predicts  a  value  of  6*  juT  —  40  at  the  RANS/LES  interface;  consequently,  we  set  Tavt  =  14  (corresponding 
to  Td  ~  42). 

One  limitation  of  the  current  approach  lies  in  the  fact  that  the  error  is  based  on  the  Reynolds  shear  stress, 
which  is  not  a  coordinate-invariant  quantity;  in  complex  geometries  or,  as  in  this  case,  if  (u'v')  is  not  the  only 
non-zero  off-diagonal  component  of  the  Reynolds  stress  tensor,  the  error  definition  might  not  be  meaningful 
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Figure  13:  3D  boundary  layer,  (a)  Streamwise  and  spanwise  skin  friction  coefficients,  C/  x  and  C/i2;  ( b ) 
mean  streamwise  velocity  profiles  and  (c)  mean  spanwise  velocity  profiles  at  the  locations  indicated  by  a 

vertical  line  in  part  (a).  o  Reference  LES; .  SA  RANS; - error  based  on  (uV): - production-based 

error. 
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(a) 


(b) 


(c) 


Figure  14:  3D  boundary  layer,  (a)  Streamwise  and  spanwise  skin  friction  coefficients,  Cf%x  and  Cj  z\  (b) 
mean  (v'wf)  Reynolds  stress  and  (c)  mean  (u'v')  Reynolds  stress  at  the  locations  indicated  by  a  vertical  line 
in  part  (a).  o  Reference  LES; SA  RANS; error  based  on  (u'v')  - prod  net  ion- based  error. 
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Figure  15:  Contours  of  streamwise  velocity  fluctuations  in  a  xz— plane  near  the  wall,  y/S*/o  —  0.18.  (a) 
Reference  LES  (only  part  of  the  domain  is  shown);  (b)  hybrid  RANS/LES. 


The  proposed  RANS/LES  merging  approach  is  predicated  on  the  RANS  calculation  being  accurate  in  the 
interface  region,  which  in  practice  restricts  the  interface  to  lie  in  a  thin,  attached  shear  layer,  where  the 
definition  of  directions  along  and  normal  to  the  shear  layer  is  unique.  Despite  this  practical  observation,  it 
would  be  desirable  to  develop  a  more  universal  error  definition,  and  in  particular  one  that  is  invariant  to  the 
choice  of  the  frame  of  reference.  One  possible  choice  to  satisfy  this  requirement  is  to  base  the  error  on  the 
production  of  turbulent  kinetic  energy,  which  satisfies  the  desired  invariance  properties.  Thus,  we  define  a 
production-based  error 

ep(y,  z,  t)  =  -{u'iv!j),ar9el (Sij)tar3et  +  <Sy)'  (22) 

(where  the  dependence  on  x0,  t/,  z  and  /  has  been  omitted,  and  StJ  is  the  strain-rate  tensor).  We  compared 
simulations  with  the  product  ion- based  error  defined  above  in  addition  to  those  with  the  standard  definition 
of  c. 

Figure  13  shows  the  skin  friction  coefficient  and  mean  velocity  profiles  in  the  streamwise  and  spamvise 
directions.  The  development  of  CfyX  is  similar  to  that  of  the  ZPG  case  (with  an  init  ial  decrease  followed  by 
a  quick  recover).  For  the  lateral  flow  we  observe  remarkably  good  agreement  with  the  reference  simulation 
beginning  from  the  end  of  the  controlled  region,  approximately  two  boundary  layer  thicknesses  downstream  of 
the  inflow.  The  principal  shear  stress,  (uV)  and  the  secondary  ones  (vfwf)  develop  quite  rapidly  (Figure  14). 
Notice  that  by  adding  the  forcing  into  the  v  equation  we  amplify  (v9v')dW/dy,  i.e.,  the  production  of  the 
secondary  stress,  as  well  as  the  production  of 

Figure  15  shows  streamwise  velocity  fluctuations  in  a  plane  parallel  to  the  wall.  We  observe  again  a  very 
rapid  build-up  of  the  streaky  structure  starting  from  the  synthetic  turbulence. 

3.3  Conclusions 

We  have  further  developed  the  controlled- forcing  method  for  turbulence  generation  at  RANS/LES  interfaces. 
The  method  has  shown  itself  to  be  robust  and  efficient.  By  studying  the  zero- pressure-gradient  boundary 
layer  (which,  despite  its  simplicity  is  one  of  the  most  challenging  cases  examined)  we  developed  guidelines 
for  the  determination  of  the  parameters  of  this  method. 

We  found  that  it  is  extremely  important  that  the  data  compared  with  the  target  Reynolds  stress  (or 
production)  is  averaged  over  a  time  interval  consistent  with  the  physics  of  the  flow.  Averaging  over  a  time 
interval  of  the  order  of  the  local  integral  time-scale  resulted  in  shorter  transient  and  more  rapid  development 
of  physically  realistic  turbulent  eddies.  The  inclusion  of  non-local  information  (as  is  the  case  when  spamvise 
averaging  is  also  applied)  is  generally  not  beneficial. 
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The  controller  constants  Kp  and  Ki  play  a  less  significant  role.  Giving  excessive  weight  to  the  integral 
error,  or  insufficient  one  to  the  proportional  part,  results  in  instabilities  of  the  flow.  For  a  wide  range  of 
values  of  Kp  and  Ki,  on  the  other  hand,  the  controller  was  stable  and  the  How  statistics  were  insensitive  to 
the  constants. 

Of  the  Hows  examined,  we  found  the  zero- pressure-gradient  and  favorable- press  lire- gradient  boundary 
layers  to  be  the  most  challenging:  in  the  zero- pressure-gradient  case  no  external  mechanism  amplify  the 
instability  of  the  How,  so  that  the  generation  of  turbulence  is  entirely  due  to  the  synthetic  turbulence 
and  controlled  forcing:  since  the  latter  is  less  effective  in  the  outer  layer  (where  the  turbulence  production 
mechanisms  are  weaker),  we  observe  the  slow  development  of  the  outer  layer.  The  favorable-pressure-gradient 
boundary  layer  presents  a  different  set  of  problems:  The  acceleration  results  in  relaminarization  of  the  How. 
Therefore,  the  injection  of  energy  due  to  the  forcing  drives  the  How  in  the  wrong  direction.  In  addition, 
RANS  models  are  not  very  accurate  in  the  strong  acceleration  region,  which  results  in  three  sources  of  error: 
Hrst,  the  mean  velocity  at  the  interface  is  far  from  the  result  obtained  with  the  full-domain  LES.  Second, 
the  target  stresses  are  not  close  to  the  “correct”  ones;  thus,  the  controller  drives  the  solution  towards  an 
incorrect  one.  Finally,  the  time-scale  obtained  from  the  RANS  is  substantially  different  from  that  obtained 
from  the  reference  calculation,  so  that  the  moving  average  is  suboptnnal. 

The  method  worked  well  in  t  he  adverse-pressuro-gradient,  case,  in  which  the  deceleration  (and,  to  an  even 
greater  extent,  the  separation)  amplifies  the  disturbances,  thus  aiding  the  generation  of  turbulent  eddies. 
This  result  is  consistent  with  the  findings  of  Terracol  [40],  who  had  some  success  using  synthetic  turbulence 
generation  in  the  ad  verse-pressure-grad  ient  region  on  the  suction  side  of  a  turbine  blade.  In  the  3D  boundary 
layer  the  turning  of  the  streamline  imposed  by  the  spanwise  pressure  gradient  drives  the  lateral  How,  with 
secondary  stresses  that  develop  from  nearly  zero  at  the  RANS/LES  interface.  The  controller  proves  to  be 
able  to  seed  the  turbulence  enough  that  the  growth  of  the  secondary  stresses,  and  hence  the  lateral  mean 
velocity  profile,  can  be  predicted  accurately. 

In  summary,  we  have  shown  that  the  controlled  forcing  method  is  robust  and  efficient ,  The  least  favorable 
results  are  obtained  in  mildly  unstable  Hows  like  the  ZPG  pressure  gradient,  or  in  re-laminarizing  ones. 
Even  in  these  situations,  within  5  boundary  layer  thicknesses  of  the  end  of  the  control  region  a  realistic  How 
was  established.  In  realistic  applications,  especially  in  curved  Hows  and  in  adverse  pressure  gradients  (an 
important  class  of  applications)  the  controlled  forcing  gives  very  rapid  development  of  realistic  eddies.  The 
application  to  a  realistic  complex  How  (the  How  over  a  turbine  blade  or  airfoil,  for  instance),  is  the  next  step 
in  the  development  of  this  method.  An  important  extension,  that  also  needs  to  be  investigated,  is  that  to 
compressible  Hows.  In  particular,  it  needs  to  be  shown  that  the  energy  added  to  the  system  through  the 
forcing  remains  as  kinetic  energy  and  is  not  transferred  to  internal  energy. 


4  Application  to  separated  flows 

4.1  Introduction 

The  target  application  of  this  research  is  the  application  of  hybrid  methods  to  Hows  (such  as  those  sketched 
in  Figure  2)  in  which  separation  may  be  unsteady  or  mild;  in  these  cases  one  would  want  to  switch  from  the 
RANS  to  the  LES  method  upstream  of  the  separation,  to  capture  the  eddies  generated  in  that  region.  In  this 
type  of  flows,  however,  an  accurate  prediction  of  the  turbulent  momentum  transport  is  extremely  important, 
since  errors  in  the  mean  velocity  profile  can  cause  premature  separation  [24,  23]).  This  phenomenon  is 
one  of  the  acknowledged  weaknesses  of  the  Detached  Eddy  Simulation  approach(see  the  recent  review  by 
Spalart  [37]).  The  rapid  generation  of  turbulent  eddies  is,  therefore,  critical. 

We  will  consider  the  flow  configuration  shown  in  Figure  16.  The  geometry  consists  of  a  Hat  plate  followed 
by  a  smoothly  contoured  ramp  and  another  Hat  plate  region.  Experiments  on  this  geometry  were  conducted 
by  Song  and  Eaton  [33].  As  the  turbulent  boundary  layer  goes  over  the  ramp,  the  flow  expands,  creating  an 
adverse  pressure  gradient  that  causes  the  How  to  separate  on  the  ramp.  The  How  subsequently  re-attaches 
on  the  flat  plate  region.  Experimental  data  are  available  at  various  streamwise  locations:  upstream  of  the 
ramp,  at  the  point  where  the  ramp  begins,  at  the  separation  point,  at  the  trailing  edge  of  the  ramp,  at 
the  re-attachment  point  and  at  two  locations  in  the  recovery  region.  Velocities  are  normalized  by  the  free¬ 
st  ream  velocity  at  the  reference  location  xp  =  —2.  Although  geometrically  simple,  the  contoured  ramp  is 
a  challenging  test  case.  Downstream  of  the  equilibrium  boundary  layer  the  adverse  pressure  gradient  and 
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Interface  region 


Figure  16:  Sketch  of  the  hybrid  RANS/LES  of  the  How  over  a  contoured  ramp. 


streamline  curvature  affect  turbulence  significantly.  The  separation  is  due  to  the  pressure  gradient,  and  was 
found  to  be  highly  unsteady  and  three-dimensional.  Incorrect  prediction  of  the  separation  point,  moreover, 
affects  the  downstream  development  of  the  flow  quite  significantly  through  the  instability  of  the  separated 
shear  layer.  As  mentioned  above,  premature  separation  is  often  the  result  of  the  application  of  hybrid 
RANS/LES  methods  to  Hows  with  shallow  separation  [37]. 

The  hybrid  calculation  uses  the  Spalart-Allmaras  [34]  turbulence  model  in  the  RANS  region.  The  same 
model  is  used  in  the  LES  region,  but  the  length  scale  is  modified  (see  below)  to  yield  an  SGS  model  of  the 
type  used  in  most  I)ES  applications.  The  interface  between  the  RANS  and  LES  regions,  shown  in  light  blue 
in  Figure  16,  is  where  the  forcing  is  applied. 

In  the  following  we  will  first  describe  the  geometry  and  the  configuration  parameters.  Then,  we  will 
present  the  equations  of  motion,  with  the  addition  of  the  forcing  terms;  we  will  then  present  the  turbulence 
models  used  for  the  calculation.  Finally,  we  will  present  and  discuss  the  results  of  the  hybrid  calculations. 


4.2  Problem  formulation 

4.2.1  Governing  equations 


The  governing  equations  are  the  Hltered  Navier-Stokes  equations,  with  the  addition  of  a  forcing  term,  /,, 
active  only  in  the  forcing  region: 


dUi 

Oxj 


&Ui 

~dt 


dUjUi 

dxj 


0 


1  d2Ui  1  dp  drtj 

Re  Oxjdxj  p  dxi  dx3  ^  1 


(23) 

(24) 


The  forcing  term  is  set  to  zero  everywhere,  except  in  a  region  immediately  downstream  of  the  RANS/LES 
interface.  It  is  designed  to  generate  velocity  fluctuations  analogous  to  those  introduced  through  the  synthetic 
turbulence  generation  described  in  Section  2.2.3,  and  to  enforce  the  controlled-forcing  technique  described 
there  Therefore,  we  let  the  force  be  given  by 

f,  =  -  ——Jjf  —  +-  <W,CF.  xif  <  X  <  +  aSref,  (25) 

where  u\  are  the  synthetic  turbulence  fluctuations  defined  by  (8-13),  ff  F  is  the  controlled  forcing  added  to 
the  ^-momentum  equation  only  and  defined  by  (15-17),  xif  is  the  location  of  RANS/LES  interface,  a  is  a 
parameter  which  defines  how  far  the  synthetic  turbulence  forcing  is  applied  from  the  RANS/LES  interface, 
and  (5ref  is  the  boundary  layer  thickness  at  the  reference  location. 
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Figure  17:  Computational  grid  used  for  the  hybrid  calculations.  Only  every  5th  grid  line  is  shown. 


4.2.2  Turbulence  model 

In  this  calculation,  the  Spalart-Allmaras  model,  described  in  Section  2.2.2,  is  used  both  in  the  RANS  and 
in  the  LES  regions.  In  the  RANS  region,  the  length-scale  d  is  given  by  the  distance  from  the  wall,  yw\  in 
the  LES  region  it  is  equal  to 

ed  —  Cdes^  (26) 

where  Goes  =  0.65  and  A  =  max(Ax,  At/,  Az).  The  resulting  model  relaxes  to  the  Smagorinsky  model  [31, 
14]  when  the  small  scales  are  in  equilibrium,  but  is  better  able  to  include  non-equilibrium  effects  through 
the  use  of  the  transport  equation  for  the  eddy  viscosity.  The  transition  between  RANS  and  LES  mode  is 
achieved  by  decreasing  the  length  scale  from  its  RANS  value  to  LES  value  linearly  over  a  streamwise  distance 
of  Src f  after  the  nominal  interface,  x  —  —2. 

4.2.3  Geometry  and  problem  parameters 

The  governing  differential  equations  (23),  (24)  and  (4)  are  discretized  on  a  noil-staggered  grid  using  a 
curvilinear  finite  volume  code.  The  method  of  Rhie  and  Chow  [25]  is  used  to  avoid  pressure  oscillations. 
Both  convective  and  diffusive  fluxes  are  approximated  by  second -order  central  differences.  A  second-order 
semi-implicit  fractional-step  [10]  procedure  is  used  for  the  temporal  discretization.  The  Crank-Nicolson 
scheme  is  used  for  the  temporal  discretization  of  wall-normal  diffusive  terms,  and  the  Adams-Bashforth 
scheme  is  used  for  the  temporal  discretization  of  all  the  other  terms.  Fourier  transforms  are  used  to  reduce 
the  three-dimensional  Poisson  equation  into  a  series  of  two-dimensional  Helmholtz  equations  in  wavenumber 
space,  which  are  then  solved  iteratively  using  the  Biconjugate  Gradient  Stabilized  (BCGSTAB)  method. 
The  code  is  parallelized  using  the  MP1  message- passing  library  and  the  domain-decomposition  technique, 
and  has  been  extensively  tested  by  Silva  Lopes  and  Palma  [29]  in  isotropic  turbulence  and  by  Silva  Lopes  ct 
al.  [30]  in  an  S-shaped  duct. 

The  geometry  used  for  the  present  calculation  is  shown  in  Figures  16  and  17.  In  the  numerical  calculations, 
all  the  lengths  are  normalized  by  the  ramp  length,  Lr  =  70  min.  The  Hat  plate  section  preceding  the  ramp  is 
2 Lr  long;  the  radius  of  curvature  of  the  ramp  is  1.814 Lr,  its  height  is  0.3 Lr.  In  the  numerical  calculations, 
the  flat  plate  region  following  the  ramp  has  a  length  of  6 Lr  followed  by  a  buffer  zone  for  the  outflow 
boundary.  The  computational  domain  height  is  1.8714 Lr  which  is  same  as  the  wind  tunnel  height.  To  save 
computational  cost,  the  upper  boundary  layer  near  the  top  wall  is  not  resolved  but  a  free  slip  boundary 
condition  is  used. 

The  grid  used  for  the  present  calculations,  shown  in  Figure  17,  uses  834x90x180  grid  points.  In  terms 
of  the  friction  velocity  at  the  reference  location  x  =  —2,  the  grid  spacing  along  streamwise  and  spanwise 
directions  are  Ax+  =  20  and  Az+  =  8.33  which  is  sufficient  for  a  wall- resolved  LES  calculation.  The  mesh  is 
refined  in  the  separation  region  and  near  the  wall.  The  Reynolds  number  of  the  calculation  was  1140  (based 
on  Ur,  the  freestream  velocity  at  the  reference  location  x  =  —  2  and  the  momentum  thickness  0). 

4.3  Results  and  discussion 

We  performed  four  calculations:  first  a  full  LES  of  the  entire  domain.  In  this  case  we  used  planes  of  data  from 
a  separate  calculation  as  inflow  condition.  Then,  we  performed  three  hybrid  cases.  In  the  first,  no  forcing 
was  applied:  ft  =  0  in  (24).  The  perturbations  present  in  the  flow  from  the  initial  condition  (which  was  one 
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Figure  18:  Flow  visualization  for  the  reference  LES  and  hybrid  calculations.  The  side  plane  shows  contours 
of  the  time-averaged  streamwise  velocity  U ,  while  the  isosurfaces  of  Q  =  60 Ur/L\  (coloured  by  the  vort.icit\ 
magnitude  |u;|)  highlight  the  vortical  structures. 


of  the  fields  of  the  full  LES)  were  allowed  to  be  amplified  and  transition  in  the  LES  region.  Another  hybrid 
case  used  the  synthetic  turbulence  but  no  controlled  forcing:  f(  F  =0  in  (25).  In  this  case  only  random 
fluctuations  were  introduced  by  the  forcing.  The  force  /,  was  active  for  a  length  Lf  —  5^rcf  downstream  of 
the  interface  (shorter  distances  were  tried,  but  were  found  to  be  insufficient  to  generate  sustained  turbulent 
eddies).  Finally,  we  performed  a  calculation  in  which  fi  included  both  synthetic  turbulence  and  controlled 
forcing.  The  forcing  in  this  case  was  active  only  for  a  length  Lf  =  3£rcf. 

Figure  18  shows  visualizations  of  the  flow  for  the  four  cases.  The  turbulent  eddies  are  visualized  bv 
isosurfaces  of  the  second  invariant  of  the  velocity  gradient. 


(where  ilij  is  the  large-scale  rotation-rate  tensor),  coloured  by  the  vorticity  magnitude. 

In  the  reference  LES  one  can  observe  that  the  How  upstream  of  the  ramp  contains  coherent  structures, 
highlighted  by  the  contours  of  streamwise  velocity  and  the  vortical  structures.  Turbulence  is  enhanced  first 
by  the  adverse  pressure  gradient,  and  further  by  the  separation,  leading  to  more  (and  stronger)  coherent 
eddies  in  the  initial  part  of  the  recovery  region.  In  the  hybrid  calculation  with  no  forcing  the  How  is  smooth 
in  the  RANS  region  (as  expected)  but  no  eddies  are  generated  in  the  LES  region  until  separation  occurs. 
The  absolute  instability  of  the  separated  shear  layer  naturally  amplifies  any  perturbation  present  in  the  flow, 
so  that  sustained  turbulence  can  be  observed  in  the  recovery  region.  The  absence  of  eddies  in  the  interface 
region,  however,  leads  to  quantitative  inaccuracies,  as  will  be  shown  below.  One  effect  of  the  forcing  is 
to  create  an  unphysical  How  in  the  interface  region,  with  excessive  generation  of  coherent  eddies  with  very 
high  vorticity  (significantly  higher  than  in  the  rest  of  the  How).  As  mentioned  above,  the  forcing  based 
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Figure  19:  Contours  of  streamwise  velocity  and  streamlines,  (a)  Full  LES:  ( b )  Hybrid  RANS/LES,  no 
forcing  at  the  interface;  (c)  Hybrid  RANS/LES,  synthetic  turbulence  forcing  at  the  interface;  (d)  Hybrid 
RANS/LES,  synthetic  turbulence  and  controlled  forcing  at  the  interface. 


on  synthetic  turbulence  needed  to  be  applied  in  a  longer  region  to  generate  sustained  turbulence.  When 
controlled  forcing  is  applied,  the  extent  of  the  linphysical  region  is  much  reduced.  For  both  cases,  however, 
the  How  structure  resembles  that  of  the  reference  LES  by  the  beginning  of  the  ramp,  and  the  separation  and 
recovery  regions  appear  realistic. 

The  effect  of  the  lack  of  coherent  eddies  that  cause  the  momentum  transport  is  evidenced  by  the  mean 
streamlines  of  the  flow,  Figure  19.  Compared  with  the  full  LES,  the  case  in  which  no  forcing  is  applied  at 
the  interface  has  early  separation  (at  xs  ~  0.22,  rather  than  xs  =  0.36,  as  predicted  by  the  full  simulation). 
The  recirculation  region  is  also  much  longer,  with  a  complex  structure  including  one  principal  and  two 
secondary  bubbles.  When  forcing  is  applied  at  the  RANS/LES  interface  the  shape  of  the  separated-flow 
region  predicted  by  the  hybrid  calculations  is  in  much  better  agreement  with  the  full  simulation.  When  the 
forcing  based  on  synthetic  turbulence  only  is  used  the  separation  is  slightly  delayed  ( xs  =  0.41  vs.  0.36), 
whereas  when  controlled  forcing  is  also  used  the  separation  point  is  predicted  correctly.  The  reattachment 
point  is  predicted  slightly  early  (xr  =  1.35)  when  synthetic  turbulence  forcing  only  is  used,  whereas  a  value 
in  agreement  with  the  full  LES  (xr  =  1.44)  is  predicted  when  controlled  forcing  is  used. 

The  qualitative  discussion  based  on  the  How  visualization  is  supported  by  the  quantitative  comparisons. 
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Figure  20:  Skin-friction  coefficient. 


The  skin-friction  coefficient 


(28) 


(where  Ur  is  the  average  velocity  at  the  reference  location  Xr  —  —2)  is  shown  in  Figure  20.  The  sudden 
decrease  of  the  Cj  at  the  RANS/LES  interface  in  the  calculation  with  no  forcing  is  due  to  the  fact  that 
changing  the  length-scale  in  the  Spalart.-Allmaras  model  decreases  the  eddy  viscosity.  Since  no  eddies  are 
present  in  this  region  (Figure  18),  the  turbulent  wall-normal  momentum  transfer  is  effectively  decreased, 
leading  to  a  more  laminar-like  profile  with  decreased  wall  stress.  The  lower  momentum  of  the  fluid  near  the 
wall  results  in  an  early  separation,  and  in  a  much  more  extended  separation  bubble.  When  the  synthetic 
turbulence  is  used,  we  observe  a  significant  increase  of  the  skin  friction  in  the  region  where  the  equations 
are  manipulated  (which  is,  of  course,  expected  to  be  unphysical),  but  a  recovery  to  the  reference  LES  values 
shortly  after  the  interface  ends.  Best  results  are  obtained  when  the  synthetic  turbulence  is  coupled  with  the 
controlled  forcing:  the  extent  of  the  unphysical  interface  region  was  reduced,  and  very  good  agreement  was 
achieved  with  the  reference  LES. 

The  profiles  of  the  mean  streamwise  velocity  and  Reynolds  shear  stress  are  shown  in  Figure  21.  When  no 
forcing  is  used,  the  lack  of  momentum-supporting  eddies  discussed  above  results  in  vanishing  resolved  shear 
stress,  which  are  not  compensated  by  the  modeled  ones,  Figure  21(b).  The  decreased  momentum  transport 
results  in  a  more  laminar-like  velocity  profile,  which  separates  early.  We  have  already  commented  on  the 
effect  of  the  early  separation  and  the  decreased  mixing  on  the  shape  of  the  separation  region  when  no  forcing 
is  used.  Note  that  even  at  x  =  2  the  calculations  with  no  forcing  at  the  interface  do  not  agree  with  the 
full  simulations  (the  simulation  predicts  that  the  flow  is  still  separated  there),  and  only  around  x  —  4,  a 
considerable  distance  downstream  of  the  ramp  the  mean  velocity  and  Reynolds  stresses  agree  with  the  full 
LES  (the  boundary  layer  thickness  upstream  of  the  ramp  is  Svv f  =  0.33). 


4.4  Conclusions 

This  work  showed  the  feasibility  and  utility  of  the  use  of  forcing  at  the  RANS/LES  interface  in  actual 
single-block  hybrid  calculations  of  flows  with  shallow  separation.  When  no  special  treatment  is  used  at  the 
RANS/LES  interface,  the  behavior  of  the  How  after  the  RANS/LES  interface  is  incorrect.  Lack  of  resolved 
eddies  results  in  a  decrease  of  the  skin  friction,  which  indicates  a  less  full  velocity  profile  near  the  wall, 
leading  to  early  separation  and  to  a  recirculation  bubble  that  is  excessively  long  and  high.  Only  a  significant 
distance  downstream  (approximately  15  reference  boundary  layer  thicknesses)  the  How  relaxes  to  its  correct 
state. 

When  forcing  is  introduced  immediately  downstream  of  the  RANS/LES  interface,  on  the  other  hand, 
momentum-transporting  eddies  are  generated  more  effectively,  the  prediction  of  the  velocity  profile  is  more 
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Figure  21:  (a)  Velocity  profiles;  (b)  Reynolds  shear  stress  profiles. 
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accurate,  and  the  separation  point  (and,  hence,  the  behavior  of  the  How  downstream)  is  predicted  well. 

The  controlled  forcing  was  particularly  effective  in  reducing  the  extent  of  the  region  in  which  the  Navier- 
Stokes  equations  had  to  be  manipulated.  The  comparison  with  the  full  LES  and  with  the  experimental  data 
was  quite  good. 

5  Summary 

We  have  demonstrated  the  applicability  of  synthetic  turbulence  coupled  with  the  controlled-forcing  method 
to  generate  eddies  at  the  interface  between  RANS  and  LES  in  hybrid  calculations.  By  investigating  the 
effects  of  the  model  parameter,  and  developing  guidelines  to  set  their  values,  we  obtained  shorter  transition 
regions  and  improved  model  accuracy.  It  was  found  that  realistic  turbulence  (with  the  correct  statistics) 
can  be  generated  within  5  boundary- layer  thicknesses  of  the  RANS/LES  interface  even  in  cases  (such  as  the 
accelerating  boundary  layer)  in  which  the  data  supplied  by  the  RANS  is  inaccurate,  or  the  assumptions  on 
which  the  forcing  is  based  are  invalid. 

We  applied  this  technique  in  a  variety  of  turbulent  Hows  including  phenomena  such  as  freestream  accel¬ 
eration,  separation,  and  mean  How  three-dimensionality,  with  uniformly  good  results.  A  single-block  hybrid 
calculation  was  then  performed  in  a  geometry  that,  although  simple,  presented  several  difficulties:  because  of 
the  shallow  pressure-driven  separation,  incorrect  prediction  of  the  upstream  How  results  in  significant,  errors 
that  propagate  downstream,  as  the  shear  layer  instability  acts  as  an  error  amplifier.  We  demonstrated  that 
unless  turbulent  eddies  are  artificially  generated  at  the  RANS/LES  interface,  very  significant  errors  appear 
in  the  How  statistics  even  at  low  order  (skin-friction  coefficient,  mean  velocity  profile).  Establishing  realistic 
turbulent  eddies  capable  to  transport  momentum  energy  and  mass,  appears  to  be  a  critical  factor  for  the 
accurate  prediction  of  shallow  separation  by  hybrid  RANS/LES  methods. 

Future  work  should  concentrate  on  the  application  of  this  method  in  compressible  boundary  layers.  In 
such  cases  the  energy  introduced  by  the  forcing  may  be  turned  into  internal  energy,  and  be  less  effective  at 
generating  turbulent  eddies,  compared  with  incompressible  cases.  The  generation  of  spurious  pressure  waves 
could  also  become  an  issue.  This  work  will  be  proposed  for  a  follow-up  grant. 
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